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We report on the results of an ongoing calculation of massive scalar one-loop one-, two-, three- and four-point 
integrals up to 0(e 2 ) which are needed in the NNLO calculation of heavy hadron production. 



1. Introduction 

The full next-to-leading order (NLO) QCD cor- 
rections to hadroproduction of heavy flavors have 
been completed as early as 1988 |1I2) . They have 
raised the leading order (LO) estimates |3] but 
were still below the experimental results (see e.g. 

In a recent analysis theory moved closer to 
experiment 0]. A large uncertainty in the NLO 
calculation results from the freedom in the choice 
of the renormalization and factorization scales. 
The dependence on the factorization and renor- 
malization scales is expected to be greatly re- 
duced at next-to- next-to-leading order (NNLO). 
This will reduce the theoretical uncertainty. Fur- 
thermore, one may hope that there is yet bet- 
ter agreement between theory and experiment at 
NNLO. 

2. General remarks 

In Fig. ^ we show one generic diagram each 
for the four classes of gluon induced contributions 
that need to be calculated for the NNLO correc- 
tions to hadroproduction of heavy flavors. They 
involve the two- loop contribution (Fig. QJi), the 
loop-by- loop contribution (Fig.[TJ)), the one- loop 
gluon emission contribution (Fig.[TJ:) and, finally, 
the two gluon emission contribution (Fig.^Ji). We 
mention that there is an interesting subclass of 
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the diagrams in Fig.^; where the outgoing gluon 
is attached directly to the loop. One then has a 
five-point function which, when folded with the 
corresponding tree graph contribution, has to be 
calculated up to 0(e 2 ). 

In our work we have concentrated on the loop- 
by-loop contributions exemplified by Fig. 
Specifically, working in the framework of dimen- 
sional regularization, we are in the process of cal- 
culating 0(e 2 ) results for all scalar massive one- 
loop one-, two-, three- and four-point integrals 
that are needed in the calculation of hadronic 
heavy flavour production. It should be clear from 
the flavour flow of the diagrams that, apart from 
the two mass scales set by the kinematics of the 
process there is only one explicit heavy mass scale 
in the problem. The integration is generally done 
by writing down the Feynman parameter repre- 
sentation for the corresponding integrals, inte- 
grating over Feynman parameters up to the last 
remaining integral, expanding the integrand of 
the last remaining parametric integral in terms 
of the dimensional parameter e and doing the 
last parametric integration on the coefficients of 
the expansion. Because the one-loop integrals ex- 
hibit infrared (IR)/collincar (M) singularities up 
to (D(e~ 2 ) one needs to know the one-loop inte- 
grals up to 0(e 2 ) since the one- loop contributions 
appear in product form in the loop-by-loop con- 
tributions. 

The aim of our project is thus to compute all 
one-loop contributions to the two processes 

1. g+g^Q+Q 
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Figure 1. Exemplary gluon fusion diagrams for 
the NNLO calculation of heavy hadron produc- 
tion 



2. q + q^Q + Q 

up to 0(e 2 ). 

Regarding the four classes of diagrams in Fig.l 
one might then say in a very loose sense that we 
are aiming to calculate one-fourth of the NNLO 
partonic contributions to heavy hadron produc- 
tion. Nevertheless, calculating this one-fourth of 
the full problem already allows one to obtain a 
glimpse of the complexity that is waiting for us 
in the full NNLO calculation. This complexity 
does in fact reveal itself in terms of a very rich 
polylogarithmic structure of the Laurent series 
expansion of the scalar one-loop integrals as will 
be discussed later on. 

In dimensional regularization there are three 
different sources that can contribute positive em- 
powers to the Laurent series of the one-loop am- 
plitudes. These are 



1 . Laurent series expansion of scalar one-loop 
integrals 

2. evaluation of the spin algebra of the loop 
amplitudes bringing in the n-dimensional 
metric contraction g^g^ = n = 4 — 2e 

3. Passarino-Veltman decomposition of tensor 
integrals involving again the metric contrac- 
tion g^g^ =n = A-2e 

Concerning the first item the 0(e 2 ) calculation 
of the necessary one-, two- and three-point one- 
loop integrals for the loop-by-loop part of NNLO 
QCD calculation have now been completed by us. 
Two of the three massive one-loop four-point in- 
tegrals have also been done leaving us with one 
remaining four-point integral which is presently 
being worked out. We hope to be able to present 
complete results on this part of the NNLO calcu- 
lation in the near future 

Concerning the last two items (spin algebra 
and Passarino-Veltman decomposition) there ex- 
ist some partial results on this part of the NNLO 
calculation which will be given at the end of this 
presentation. 

Apart from the present discussion of NNLO 
contributions to heavy hadron production the cal- 
culation of massive scalar loop integrals up to a 
given positive power of e is of interest also in 
other contexts. For example, if the one-loop inte- 
grals appear as subdiagrams in a given divergent 
Feynman diagram one again needs to avail of the 
positive e-powers of the subdiagram. This is of 
relevance for the calculation of two-loop counter 
terms. Another example is the reduction of a 
given set of loop integrals to master integrals by 
the integration-by-parts technique jHj ■ In the re- 
duction one may encounter explicit inverse pow- 
ers of e which implies that one has to evaluate the 
master integrals up to positive e- powers. 

3. Laurent series expansion of scalar one- 
loop integrals 

In Table 1 we provide a list of the one-loop 
scalar one-, two-, three- and four-point integrals 
that need to be evaluated up to 0(e 2 ) in NNLO 
heavy hadron production. In column 2 the in- 
tegrals are identified using the notation of 
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Table 1 

List of scalar n-point functions up to C(e 2 ) 





Nomenclature of Beenakker et al. |2] 


Our nomenclature 


Novelty 


Comments 


1-point 


A(m) 






Re sj 


2-point 


B(p4 — P2,0, m) 


Si 




Re yj 




B(p3 +p4,,m,m) 


Sa 




Re, Im yj 




B(p4,0,m) 


Bs 




Re V 




B(p2, m, m) 






Re yj 




J B( P3 +P4,0,0) 


Bs 




Re, Im y 


3-point 


C(p 4 ,P3,0, m,0) 


Ci 


new 


Re, Im yj 




C(p4, -p2,0,m,m) 


C 2 


new 


Re V 




C(-p 2 ,P4,0,0,m) 






Re yj 




C(-p 2 , -pi, 0,0,0) 


c; 4 




Re, Im y 7 




C(-p2,—pi,m,m,m) 


c 5 


new 


Re, Im yj 




C(p 3 ,p 4 ,m,0,m) 


c 6 




Re, Im y 


4-point 


D(p4, -pa, -pi , 0, m, m, m) 


Di 


new 


Re, Im 




O(-p 2 ,P4,P3,0,0,m,0) 


D 2 


new 


Re, Im yj 




D(-p 2 ,p4, — pi, 0, 0, m, m) 


Da 


new 


Re yj 



When writing down a Laurent series expansion of 
these integrals in terms of e-powers one needs to 
introduce a short-hand notation for the integrals 
in order to keep the notation managable. Our 
short-hand notation for the integrals appears in 
column 3. In column 5 we comment on their real- 
ity property and tick off those integrals that have 
been completed up to now. As mentioned before 
all one-loop integrals have been done except for 
the four-point integral D\ involving three mas- 
sive propagators. We have compared our results 
to results in the literature whenever possible and 
when these were accessible. We have found agree- 
ment. In column 4, finally, we indicate which of 
our results are new. 

We mention that one also needs the imaginary 
parts of the amplitudes since the square of the 
amplitude contains also imaginary parts accord- 
ing to 



\A\ 2 = (ReAf + {ImAf 



(1) 



Note, though, that the imaginary parts are only 
needed up to 0{e 1 ) since the IR/M singularities 
in the imaginary parts of the one-loop contribu- 
tions are of 0(e~ 1 ) only. 

The scalar four-point integrals are the most 
difficult to calculate. They contain a very rich 
structure in terms of polylogarithmic functions. 
For example, the e 2 -coefficients of the Laurent 



series expansion of the four-point integrals con- 
tain logarithms and classical polylogarithms up 
to order four (i.e. Li^) (in conjunction with the 
^-functions £(2,3,4)) and products thereof, and 
a new class of functions which are now termed 
multiple polylogarithms [Jj. A multiple polylog- 
arithm is represented by 



Limit j- ■ ip^h -i ) 



dt 



dt 



dt 
o I — o 



X2X3. ..Xk — t \ t 

dt 





m 2 — 1 



O ... O — O 

X 3 ...Xk -t V t 



dt V 1 ^ 1 dt 



l-t ' 

where the iterated integrals are defined by 



dt dt 
o ... o 



a n - 1 



a-i — t 



A t 
dt n 



dt,, 



1 2 



0*71 J ^n— 1 in— 1 





dt 



Instead of using the multiple polylogarithms of 
Goncharov we have chosen to write our results in 
terms of one-dimensional integral representations 
given by the integrals 
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nl d ln(ai + cry/) h\(q 2 + cr 2 y) ln(a 3 + cr 3 y) 
) « 4 + y 



and 



ln(ai + criy)Li 2 (a 2 + a 3 y) 



where the <7i take values ±1 and the Oj's are 
combinations of the kinematical variables of the 
process. The numerical evaluation of these one- 
dimensional integral representations are quite sta- 
ble. 

A comment on the length of our expressions is 
appropiate. The untreated computer output of 
the integrations is generally quite lengthy. The 
hard work is to simplify these expressions. We 
have written semi-automatic computer codes that 
achieve the simplifications using known identities 
among poly logarithms. All intermediate steps 
have been checked numerically. Even after sim- 
plification, the final expressions are generally too 
long to be given in this presentation. As men- 
tioned earlier on, we have compared our results 
to the results of other authors when they were 
available and have found agreement. 

4. Laurent series expansion of one loop 
amplitudes 

In order to obtain the Laurent series expansion 
of the full one-loop amplitude one has to combine 
the Laurent series expansion of the scalar one- 
loop integrals with the ^-expansion from the spin 
algebra calculation and the Passarino-Veltman 
decomposition. This is the subject of a second 
paper which we are working on |HJ. As an exam- 
ple of such an amplitude result we present results 
on the contribution of the one-loop gluon trian- 
gle graph to the process g + g — > Q + Q. We 
shall only list the e— and e 2 -terms for the triangle 
graph contribution since the finite and divergent 
pieces have been calculated before (see e.g. |5]). 

Let us begin by defining a reduced amplitude 
which is obtained by extracting the spin wave 
functions from the full amplitude M, i.e. 



M = e( Pl )e( P2 )u(p 3 )M^u{ Pi ). 



(2) 



For the reduced gluon triangle amplitude 
M^s(g) with gluons and ghosts inside the tri- 
angle loop one obtains 



M^ig) = N c {-3s(B^[207bP + 12fl£ 



,(0) 



+54Ci 1) s 



6i(T a T b -i*T a ^> 1 [<f v 



[9BP-12Bi 0) + 9Ci 1] s- 8]/s +8pft# [3B { 5 0) 
+2]/s 2 )) - £ 2 (S^[621S^ 3 + 36fl£ X) + 24^ 0) 



162Cf h + 16] + &i(T a T b - T 6 T a )^i 



\27B, 



(2) 



36B 5 



24B 

>(o) 



(o) 



27Cf's-lQ]/s 



+8p£#[9B^ + 6B?> + 4]/ s 2 ))}/324. (3) 

The one- loop triangle graph amplitude M^Jg) 
has been written in terms of the Laurent series ex- 
pansion of the two massless scalar two-point and 
three-point one-loop integrals B 5 and C5 listed 
in Table 1, where the Laurent series expansions 
are defined by 



B 5 = iC £ (m 2 ) 



and 



C 4 = iC e {m 2 ) 



N 



N 
n=-2 



(4) 



(5) 



where N = 2 in our application. The Born term 
amplitude B% v appearing in is defined by [5j 

B% u = 2i{T a T b ^T b T a )(g^ 1 +p^^p^)/s, (6) 

and the coefficient C £ (m 2 ) is given by 

,2\ 2 



CM 1 ) 



r(i 



(4tt) 5 



(7) 



After insertion of the appropiate coefficient 
functions B^ and cjj^ our results can be seen to 
fully agree with the results of the authors of ^H] 
who also calculated the gluonic one-loop correc- 
tions to the three-gluon vertex with one off-shell 
gluon. 

Eq.JJJ gives an impression of the interplay of 
the three different sources of positive e-powers in 
the amplitude calculation mentioned earlier on. 
Note in particular that different orders of the Lau- 
rent series coefficients of the scalar integrals enter 
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at each order of the Laurent series expansion of 
the full amplitude. 
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